Spin Dynamics for the t-J Model 
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The spin dynamics at the finite temperature for the t-J model in the underdoped and optimal 
doped regimes is studied within the fermion-spin theory. It is shown that the dynamical spin 
structure factor spectrum at the antiferromagnetic wave vector Q — are separated as low- 

and high-frequency parts, respectively, but the high-frequency part is suppressed in the dynamical 
susceptibility spectrum x (Qi^X while the low-frequency part is the temperature dependent, which 
are in qualitative agreement with the experiments and numerical simulations. 
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The t-J model provides one of the most interesting 
and simplest models to study the physics of the doped 
Mott insulator 1 . It was originally introduced as an ef- 
fective Hamiltonian of the Hubbard model in the strong 
coupling regime 2 , where the on-site Coulomb repulsion 
U is very large as compared with the electron hop- 
ping energy t, and in this case the electrons become 
strongly correlated to avoid the double occupancy, i.e., 
Xio-C^C". < 1. The interest in the two-dimensional 
(2D) t-J model was stimulated by many researchers' sug- 
gestion that it contains the relevant physics of the copper 
oxide superconductors . 

During the past ten years, consistent experiments pic- 
ture of the magnetic properties of the copper oxide ma- 
terials have been emerged 3 ' 4 : the undoped copper ox- 
ide materials are antiferromagnetic Mott insulators, and 
upon doping with holes in the copper oxide sheets, the 
antiferromagnetic long-range-order (AFLRO) is replaced 
by short-range spin correlations, and the dynamic anti- 
ferromagnetic correlations persist even into the supercon- 
ducting state. A series of the neutron-scattering study 5,6 
of the copper oxide materials shows that there is an 
anomalous temperature dependence of the spin fluctu- 
ations near the antiferromagnetic zone center in the un- 
derdoped and optimal doped regimes, which has a much 
larger effect on the spin dynamics and leads to the un- 
usual temperature dependence of the spin-lattice relax- 
ation time. On the theoretical side, the spin dynamics 
in the underdoped and optimal doped regimes has been 
studied within the phenomenological model of antifer- 
romagnetic correlated spins 7 , and the framework of the 
phenomenological marginal Fermi- liquid theory 8 . More- 
over, the most reliable results of the spin dynamics have 
been obtained by the numerical simulations 9 and high- 
temperature series expansion 10 based on some strongly 
correlated models. It is believed that the role played by 
the magnetism, particularly the nature of spin fluctua- 
tions, is the central issue of the copper oxide materials. 
In this paper, we employ the t-J model to study the 
spin dynamics of the copper oxide materials in the un- 
derdoped and optimal doped regimes within the fcrmion- 
spin theory 11,12 , and show that our theoretical results are 
qualitative consistent with the experiments 5 and numer- 



ical simulations 9 . 

In order to account for the real experiments under the 
t-J model, the crucial requirement is to impose the elec- 
tron on-site local constraint for a proper understanding 
of the physics of the copper oxide materials 13 . For in- 
corporating the local constraint, the fermion-spin theory 
based on the charge-spin separation has been proposed 11 . 
In this approach, the constrained electron operators are 
decoupled as = h^S^ and d± = hfS^, with the 
spinlcss fermion operator hi keeps track of the charge 
(holon), while the pseudospin operator Si keeps track 
of the spin (spinon). The advantage of the fermion- 
spin theory is that the electron on-site local constraint 
for the single occupancy is satisfied even in the mean- 
field approximation (MFA) . Within the fermion-spin the- 
ory, the mean-field theory 12 in the underdoped and opti- 
mal doped regimes has been developed, where the mean- 
field order parameters are defined as % = {^t ^i+ n ) = 

(<% Si+T])' Xz — (SiSf +v ), C — (^/Z 2 )J2r 1 ,r 1 '(^i+7j^i+r 1 ')^ 

c z = (i/^)E w^fWSJW), and * = with 

fj = ±i, ±?/, and Z is the number of nearest neighbor 
sites. This mean-field theory 12 has been applied to study 
the electron spectrum, electron dispersion, and electron 
density of state of the copper oxide materials, and the 
results are in agreement with the experiments and nu- 
merical simulations. In the fermion-spin representation, 
the t-J model is written as H — H t + Hj with 

H t = -tJ2hiht +v (S+Sr +v + Srsl n ) + h.c. 
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where J e ff = J[(l — S) 2 — </> 2 ], and /x is the chemical po- 
tential. Based on the Ioffe-Larkin combination rule 14 , 
the charge dynamics of the copper oxide materials in 
the underdoped and optimal doped regimes has been 
discussed 15 by considering the charge fluctuations around 
the mean-field solution. However, the spin fluctuations 
couple only to spinons and therefore no composition law 
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is required 14 in discussing the spin dynamics, but the 
effect of holons still is considered through the holon's or- 
der parameter <j> entering in the spinon propagator. For 
discussing the spin dynamics, we need to calculate the 
second-order correction for the spinon by going beyond 
MFA. The second-order spinon self-energy diagram from 
the holon pair bubble is shown in Fig. 1. The mean- 
field spinon Green's functions D(°\k,u) and D^\k,u>), 
and mean-field holon Green's function g(°'(k,(j) have 
been given in Ref. 12 . Since the spinon operator obey 
the Pauli algebra, it is needed to map the spinon opera- 
tor into the CP 1 fermion representation or the spinless- 
fcrmion representation in terms of the 2D Jordan- Wigner 
transformation 16 for the formal many particle perturba- 
tion expansion. After such formal expansion, the spinon 
Green's function in the spinon self-energy diagram shown 
in Fig. 1 is replaced by the mean-field spinon Green's 
function D 1 - ' (k, uS), then the second-order spinon self- 
energy is evaluated as, 



S< a) (*,w) = -(2i) 2 ^(7 t - P +l P > +P+ k?B k+p , 
pp' 

Fi(k,p,p') 
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where F 1 (k,p,p') = n F (^ p+p >)[l - n F (^ p )} + [1 + 
n B (uk+ P >)}[nF(€ p ) - «f(W)]' 

F 2 (k,p,p') = n F (£ p+ j/ ) [1 - n F (£„)] - n B {uk+ P > ) [n F (£p) - 
npij^p+pi)] n F (£k) and nsi^k) are the Fermi and Bose 
distribution functions, respectively, 7^ = (1/Z) e lk ''', 
B k = ZJ eff [{2e Xz +xhk~(ex + 2x z )]/u(k), e = 1 + 
2tcj)/J e ff, and the mean- field holon excitation spectrum 
£fc and mean-field spinon excitation spectrum ujk are 
given in Ref. 12 . In this case, the full spinon Green's 
function is obtained as D~ 1 (k,uj) = D^~ 1 {k,w) — 

Ei 2) (fc^) ; 

As manifestation of the spin dynamics, the dynamical 
spin structure factor S(k, ui) and dynamical susceptibility 
x(k,uj) are given as 
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FIG. 1. The spinon's second-order self-energy diagram. 
The solid and dashed lines correspond to the holon and spinon 
propagators, respectively. 



X "(k,u;) = (l-e-^)S(k,u ) ) 
= 2lmD(k,oj), 
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respectively. Although the detailed magnetic properties 
depend on the sample preparation method in the ex- 
periments as well as the precise value and homogene- 
ity of the oxygen content, the anomalous temperature 
dependence of the spin fluctuations near the antiferro- 
magnetic wave vector Q = (it, ir) seems common 5-7 . We 
have performed the numerical calculation for the dynam- 
ical structure factor (3) and dynamical susceptibility (4). 
The result of the S(Q, u) spectra at the doping S = 0.06 
(solid line) and 6 — 0.15 (dashed line) with the temper- 
ature T = 0.2J are plotted in Fig. 2. It is shown that 
in the underdoped and optimal doped regimes there are 
the coexistence of the low- and high-frequency fluctua- 
tions in the S(Q, ui) spectra, the excitations are remark- 
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FIG. 2. The dynamical spin structure factor spectra 
S(Q,uj) in the temperature T = 0.2 J for the doping 8 = 0.06 
(solid line) and 8 = 0.15 (dashed line). 
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FIG. 3. The dynamical spin susceptibility spectra 
X (Q,w)/w at the doping 8 — 0.15 for the temperature 
T = 0.2 J (solid line), T = 0.3 J (dashed line), and T = 0.45 J 
(dotted line). 
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ably sharp, and the spectra are changed with dopings, 
which is consistent with the experiments 5 and numerical 
simulation 9 . The low-frequency peak in Q point is due to 
the antiferromagnetic fluctuations, which will be in ex- 
istence even in the undoped case, and dominate dynam- 
ical susceptibility and the neutron-scattering processes, 
while the high-frequency peak may come from the con- 
tribution of the free-fermion-like component of the sys- 
tems, which induces the main effect to the large extent 
the static spin correlation. In correspondence with the 
S(Q,uj) spectra, the numerical results of the dynamical 
susceptibility spectra \ & t the doping 5 = 0.15 

for the temperature T = 0.2 J (solid line), T = 0.3 J 
(dashed line), T = 0.45 J (dotted line) are plotted in 
Fig. 3. Comparing with the Fig. 2, it is shown that 
although the high-frequency peak is suppressed in the 
X (Q,u)) spectra, however, it still is separated from the 
low-frequency part at the antiferromagnetic wave vec- 
tor Q. Our results also indicate that the low- frequency 
peak of the dynamical susceptibility \ (Q^) is tem- 
perature dependent, while the high-frequency part is 
almost temperature independent, which are consistent 
with the experiments 5 and numerical simulations 9 . The 
present theoretical results have been used 17 to extract 
the integrated susceptibility and spin-lattice relaxation 
time, and the results shown that the integrated suscep- 
tibility exhibits the particularly universal behavior as 
l/NJ2 k x"(k,uj) cx arctan[aiw/T + a 3 (w/T) 3 ] and the 
spin-lattice relaxation time is weakly temperature depen- 
dent. 

In summary, we have studied the spin dynamics of the 
t-J model in the underdoped and optimal doped regimes 
within the fcrmion-spin theory. It is shown that there are 
two peaks for the dynamical spin structure factor at the 
antiferromagnetic wave vector Q, but the high-frequency 
peak is suppressed in the dynamical susceptibility spec- 
tra x which are qualitative consistent with the 
experiments and numerical simulations. 
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